Vortex configurations and dynamics in elliptical pinning sites for high matching fields 



in 
o 
o 

o 

Q 

m 

o 
o 

Oh' 
o3 



i 

G 

O 

o 



> 

00 

o 
o 

o 

\o 
o 

-I— > 



C 

o 
o 



X 



C.J. Olson Reichhardt 1 , A. Libal 1 ' 2 , and C. Reichhardt 1 

1 Theoretical Division and Center for Nonlinear Studies, 
Los Alamos National Laboratory, Los Alamos, New Mexico 87545 
2 Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556 

(Dated: February 6, 2008) 

Using numerical simulations we study the configurations, dynamics, and melting properties of 
vortex lattices interacting with elliptical pinning sites at integer matching fields with as many as 
27 vortices per pin. Our pinning model is based on a recently produced experimental system [G. 
Karapetrov et ah, Phys. Rev. Lett. 95, 167002 (2005)], and the vortex configurations we obtain 
match well with experimental vortex images from the same system. We find that the strong pinning 
sites capture more than one vortex each, and that the saturation number of vortices residing in a 
pin increases with applied field due to the pressure from the surrounding vortices. At high matching 
fields, the vortices in the interstitial regions form a disordered triangular lattice. We measure the 
depinning thresholds for both the x and y directions, and find distinctive dynamical responses along 
with highly anisotropic thresholds. For melting of the vortex configurations under zero applied 
current, we find multi-step melting transitions in which the interstitial vortices melt at a much 
lower temperature than the pinned vortices. We associate this with signatures in the specific heat. 

PACS numbers: PACS numbers: 74.25. Qt 



The immobilization of superconducting vortices at pin- 
ning sites is a critical area of study for device applica- 
tions, since the dissipation that occurs when the vor- 
tices move under an applied current sets a limitation on 
the effectiveness of the superconductors for lossless trans- 
port. Any superconductor contains naturally occurring 
pinning which may take the form of nonsuperconduct- 
ing inclusions, twin boundaries, or other types of defects. 
Over the past few decades, there has been a tremendous 
amount of progress in the fabrication of artificial nano- 
engineered pinning, and particularly lattices of pins. Due 
to the possibility of commensuration effects between the 
vortex lattice and the pinning lattice, the artificial pin- 
ning can be tremendously more effective at immobilizing 
the vortices and enhancing the maximum critical current 
than the naturally occurring disorder in a sample. Sig- 
nificant effort has been focused on developing and study- 
ing different types of periodic pinning, beginning with 
one-dimensional pinning created through thickness mod- 
ulation of the superconducting film This was soon 
followed by the fabrication of two-dimensional arrays of 
holes [111 HQ. 

Large holes can confine more than one flux quantum 
each [f| Q . Multiple vortices inside a single hole coalesce 
into a single multiquanta vortex, since there is no su- 
perconducting material to support distinct cores for each 
vortex. In contrast, in blind holes, which do not pass 
completely through the superconducting material, indi- 
vidual vortex cores retain their identity within the hole. 
Multiple occupation by single quantum vortices in large 
blind holes was imaged experimentally using a Bitter dec- 
oration technique Q and later studied in simulation [t| . 

Continuing advances in nanolithography focused at- 
tention on the production of small, closely spaced holes 
that could each capture a single vortex. Combined with 
this, significant advances in imaging techniques have con- 



tributed to our understanding of the interaction between 
vortices and the periodic pinning. In particular, if the ap- 
plied magnetic field exceeds the matching field at which 
the number of vortices present equals the number of pin- 
ning sites, then in the case of small holes, the extra vor- 
tices above the matching field cannot occupy the holes 
but must instead sit at interstitial locations in the super- 
conducting material surrounding the holes. These inter- 
stitial vortices are weakly pinned via the repulsion from 
vortices trapped at pinning sites, and their configurations 
have been imaged using electron-beam micro scopy t ech- 
niques 01 > scanning Hall probe microscopy [ill [T2l IL^ ] 
and scanning /iSQUID devices [LJ|. The interstitial vor- 
tex arrangements are sensitive to the structure of the 
periodic pinning, and a wide range of different types of 
artificial pinning lattices have been tested experimentally 
|l5l HH ITtI Il8| , including both symmetric pinning arrays 
as well as asymmetric arrays such as rectangular arrange- 
ments of holes |2(| . The weak pinning of the inter- 
stitial vortices can produce asymmetric critical currents 

mm mm 

Simulations of vortices in periodic arrays of singly oc- 
cupied pins give good agreement with experimental ob- 
servations and indicate the importance of the interstitial 
vortices in determi ning the dynamical response of the 
system (2^, [2(| [2?], l28| . This is particularly important 
in understanding the field dependence of the behavior 
of nanopatterned superconductors. Vortices are strongly 
pinned at the matching field when the number of vor- 
tices equals the number of artificial pins, but the effec- 
tiveness of the pinning is dramatically reduced whenever 
defects in the form of interstitials or vacancies are intro- 
duced by shifting the applied field away from its matching 
value. The result is experimentally observable as oscil- 
lations in critical current [2^ or resistivity [3(], , 
as well as through separate depinning transitions for in- 
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FIG. 1: Vortex configurations at the first nine matching fields. Black dots indicate vortex positions and pinning sites are 
indicated by ellipses. The pinning troughs connecting the pins in the y direction are not directly indicated in the figure but 
can be identified based on the vortex positions. B / — (a) 1, (b) 2, (c) 3, (d) 4, (e) 5, (f) 6, (g) 7, (h) 8, and (i) 9. 



terstitial vortices [33l HE • Due to the fact that 
the pinning of interstitial vortices is weak, pins which 
capture multiple vortices are much more effective at in- 
creasing the critical current over a wider range of field 
[33. There is, however, a trade off between saturation 
number of the pin and size of the nanoarrays; for exam- 
ple, in an imaging experiment for relatively large dots, 
as many as 280 vortices were captured by each artificial 

pin m. 

In the interest of combining small size with greater 
pinning effectiveness, recent work has shifted away from 
simple circular holes or antidots to asymmetric pinning 
sites. Asymmetry can be introduced by placing sym- 
metric pins in asymmetric unit cells containing two pins 
each |3!J |4(| ; alternatively, the pins themselves may be 
made into asymmetric shapes, such as triangles |4l|, |43 
or rectangles . Arrangements of rectangular pins are 
of particular interest since they may be used as the basis 
for vortex logic devices, in which each elongated pin can 
hold a single bit of information [3. EBf . 



A new fabrication technique has recently been devel- 
oped to create well controlled elongated pinning sites in a 
sample with an atomically flat surface |46l |47| . This per- 
mits direct imaging of the vortex positions using scanning 
spectroscopy methods. Due to the extremely high reso- 
lution of this imaging technique, vortex configurations 
can be observed at fields well above matching, showing 
interesting effects such as the formation of a disordered 
triangular vortex lattice in the region between pinning 
sites. Obtaining an image with a scanning technique 
is a time-consuming process, so dynamical information 
about the vortices is not available from the experiments 
of Refs. 0, 0|- K would be interesting to explore the 
dynamical behavior of the vortices under applied currents 
or finite temperatures at high matching fields. Existing 
simulations have been confined to low matching fields, 
with studies reaching at most the ninth matching field 
|27|. and higher matching fields have never been consid- 
ered. The particular pinning geometry used in Ref. 0] 
also has unique aspects which have not been explored 



3 



previously. 

In this work, we study the static and dynamic behav- 
iors of vortices interacting with elongated pinning sites 
of the type created in Ref. [4(| for matching fields as 
high as 27 vortices per pinning site. We explore the dy- 
namic response of the vortices to applied currents, and 
find strongly asymmetric depinning forces. We also study 
the melting of the vortex configurations and find inter- 
esting two-stage melting behavior in which the intersti- 
tial vortices melt at much lower temperatures than the 
pinned vortices. The particular model of the pinning that 
we use has not been considered before and gives good 
agreement with the vortex configurations obtained in ex- 
periment. We also study much higher matching fields 
than have been simulated previously. The dynamic vor- 
tex responses that we observe are related to the particular 
symmetries present at each matching field. 

We model a two-dimensional system of N v vortices in a 
superconducting crystal containing N p artificial pinning 
sites. We assume periodic boundary conditions in the x 
and y directions. The overdamped equation of motion 
for an individual vortex i is 

^ = fr +ff + i? + {? (1) 

The damping constant r\ — (f^d/2ii^ 2 p^ in a crystal of 
thickness d. Here, 4>q = h/2e is the flux quantum, £ is 
the superconducting coherence length, and pN is the nor- 
mal state resistivity of the material. The vortex-vortex 
interaction force is 

N v 

fr = E/" A 'i(T) fij ' (2) 

where K± is the modified Bcssel function appropriate for 
stiff three-dimensional vortex lines. A is the London pen- 
etration depth, fo = 4>q/ (2-KpoX 3 ), and is the distance 
between vortices i and j. We measure time in units of 
t = rj/fo. For a NbSe2 crystal 0.1 mm thick, we take 
r] = 2.36 x 1CT 11 Ns/m, f = 6.78 x 10~ 5 iV/m, and 
t = 0.35 /is. The thermal force if arises from ran- 
dom Langevin kicks with the properties (if) = and 
(if \t)ij '(£')) = 2-qk B T8(t - £')<%• The conversion from 
physical temperature to simulation units is given by 

Using 5t = 0.01, a temperature of 4.2K corresponds to 
fx = 4.12 x 10~ 7 /o. Thus, for our measurements of trans- 
port characteristics, we take fr = 0. 

The quenched disorder if is modeled as composite pin- 
ning sites constructed to match the experimental geom- 
etry. The sample contains parabolic troughs 0.5A wide 
that run in the y direction and are spaced 2.5A apart in 
the x direction, corresponding to 0.5^tm in NbSe2. Each 
trough is decorated with parabolic traps that confine only 
in the y direction, which define the locations of the pin- 
ning sites shown in the experimental images. These traps 



are 1.5A long in the y direction and equal to the trough 
width in the x direction. The long troughs represent 
the fact that the crystal layer in Ref. ^(J is extremely 
thin between the pins in the y direction. This produces 
an energetically favorable location for a vortex. Within 
the pinning sites, the order parameter is completely sup- 
pressed by the gold in the actual experimental sample, 
and the vortices inside the pinning site would form a 
macrovortex state which is, however, not circular due to 
the pinning geometry In our model, macrovortex forma- 
tion is not permitted, but the configuration of vortices 
inside a given pin approximates the noncircular current 
configuration expected to occur in the experimental ge- 
ometry We use pinning strengths of 12.0 for the deco- 
rating troughs and 90.0 for the channel. This produces a 
square array of pinning sites each of which is elongated 
in the y direction, as illustrated in Fig. ^ 

The Lorentz driving force from an external applied cur- 
rent is given by i d , which may be applied either in the x 
direction, i d — /^x, or in the y direction, i d = fgy. We 
measure the corresponding velocity responses in either x, 

V x = (l/JVt.XE^iVi-x), ory, V y = (l/^)(E 4 ="i ^ -y). 
The system size is measured in units of A and the forces 
in terms of fo. Most of the results presented here are for 
a sample of size 10A x 10A containing N p = 16 pins. We 
consider filling fractions up to B/B^ — 27, corresponding 
to 432 vortices. At these fields, even though the vortex- 
vortex interaction is of a short-range form appropriate 
for a crystalline sample, many vortices fall within the 
interaction range of a given vortex, and as a result the 
simulations become computationally intense. It is for 
this reason that we consider the sample size illustrated 
in the figures. We have also tested larger samples of size 
30A x 30A containing 144 pinning sites and find substan- 
tially identical results. 

We first consider the vortex arrangements at integer 
filling fractions B/B^ = N v /N p . For each field, we ini- 
tialize the sample at a high temperature so that the vor- 
tices are diffusing freely, and then slowly decrease the 
temperature in order to anneal the system toward the 
T = ground state. Figure 2] illustrates the zero tem- 
perature vortex configurations at the first nine match- 
ing fields for a sample with N p = 16 pinning sites. For 
B/Bc/t — 1 to 4, the vortices are confined fully within the 
pinning sites and form single lines. At B/B^ = 5, the 
pinning confinement is not strong enough to capture a 
fifth vortex inside the pin, so one vortex per unit cell oc- 
cupies an interstitial position. If the pinning force existed 
only within the ellipses drawn in the figure, the intersti- 
tial vortex would occupy the low energy point at the same 
y position located halfway between rows of pins, but it 
would move over by half a pinning lattice constant in or- 
der to sit midway between columns of pinning. This does 
not happen due to the presence of the additional pinning 
trough introduced in our model. A similar vortex con- 
figuration for one interstitial vortex per unit pinning cell 
is observed in the experiment due to the nonuniform ion 
beam patterning of the pinning geometry in the x and y 
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FIG. 2: Vortex configurations at the second nine matching fields. Black dots indicate vortex positions and pinning sites are 
indicated by ellipses. B/Bj, = (a) 10, (b) 11, (c) 12, (d) 13, (e) 14, (f) 15, (g) 16, (h) 17, and (i) 18. 



directions, and the configuration in Fig.^e) matches well 
with the configuration shown in Fig. 2(b) of Ref. ^(|. At 
B j ' = 6, a second interstitial vortex can not fit inside 
the pinning trough, so it occupies the next most favor- 
able location centered between four pinning sites. The 
result is the formation of rows of interstitial vortices, and 
Fig. njf) resembles the configuration shown in Fig. 2(c) 
of Ref. 0]. At B/B^ > 6, we begin to observe some 
disorder in our vortex configurations, and for B / ' B$ = 6, 
a few of the pins have captured five vortices instead of 
four. This trend continues at higher fillings: the occu- 
pation number of the pins does not saturate but slowly 
increases as the overall higher vortex density outside the 
pinning site helps to stabilize a larger number of vortices 
inside the pinning site. A gradual increase in hole sat- 
uration number was also observed in the experiments of 
Ref. 0| . At B/B^ = 7 the interstitial vortices begin to 
spread out into the free channel between columns of pin- 
ning sites, and this continues at B/B^ = 8, where all pins 
now capture five vortices. The interstitial vortex chan- 
nels in Fig. [Jh) match well with Fig. 2(d) of Ref. 



By B/B^ = 9, some defects occur in the interstitial chan- 
nels and the vortices start to buckle out along the space 
between pinning rows in order to form a second intersti- 
tial column. Figures ^i) is in good agreement with the 
configuration obtained in Fig. 3(a) of Ref. |46j . 

Vortex configurations at the second nine matching 
fields, B/Bj, = 10 to 18, are illustrated in Fig. [3 For 
BjB^ = 10 and 11, we find a gradual buckling tran- 
sition of the columns of interstitial vortices filling the 
vertical gaps between the pins. The vortex configura- 
tion in the columns gradually changes from one vortex 
wide at B/B^ = 10 in Fig. 01a) to two vortices wide by 
B/B^ = 15 in Fig. HT). At B/B, p = 16 in Fig. Ejfg) 
we obtain a configuration that agrees well with Fig. 3(b) 
of Ref. |4^| . The buckling transition begins anew as the 
field is further increased, and by B/B^ = 18 in Fig. Eti) 
portions of the interstitial channels contain a vortex ar- 
rangement three vortices wide. The pinning sites con- 
tinue to capture larger numbers of vortices as the field 
increases, and at BjB^, = 15 and above these configura- 
tions also begin to buckle inside the pinning sites. The 
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FIG. 3: Vortex configurations at the third nine matching fields. Black dots indicate vortex positions and pinning sites are 
indicated by ellipses. B/Bj, = (a) 19, (b) 20, (c) 21, (d) 22, (e) 23, (f) 24, (g) 25, (h) 26, and (i) 27. 




FIG. 4: The critical depinning force for two different driving 
directions at matching fields B/B^ — 1 to 27. Circles: f* for 
driving in the x direction; squares: f% for driving in the y 
direction. 




FIG. 5: Representative velocity-force curves for x-direction 
driving, V x versus /£, for B/B^ — (from right to left) 1, 9, 
11, 14, 17, 20, 22, 24, and 27. 
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FIG. 6: Illustration of the vortex configurations just before 
depinning under application of a driving force in the x direc- 
tion (to the right in the figure). Black dots indicate vortex 
positions and pinning sites are indicated by ellipses; the pin- 
ning troughs between the pins running in the y direction are 
not directly indicated in the figure. B / = (a) 7, (b) 9, (c) 
11, (d) 17, (e) 20, and (f) 27. 



small portion of the pinning channel between the pins 
in the y direction also captures more than one vortex 
starting at Bj 'B$ = 13. 

We illustrate vortex configurations at matching fields 
much higher than have been simulated previously, 
through B/Btf, = 27, in Fig. [3] Vortices continue to pop- 
ulate the interstitial region of the channels as the field 
is further increased, and the vortex occupation number 
of the pins and of the small trough region connecting 
the pins continues to slowly increase with B/B^. We do 
not find perfectly parallel rows of vortices as was indi- 
cated schematically in Fig. 3(c) of Ref. ^fji but instead 
observe a marked modulation of the vortex density in- 
side the channel caused by the strong repulsion from the 
large number of vortices confined within each elongated 
pin. This is in general agreement with the actual image 
in Fig. 3(c) of Ref. ..46], which compares well with our 
Fig.[3fe) at B/B^ = 23. There is still a tendency for the 
interstitial vortices to form rows parallel with the x axis 
and passing between pinning rows, which is particularly 
pronounced at B/B^ — 20 in Fig. I^fb) and to a lesser 
degree at B/B^ = 22 in Fig.[3Jc). As we will show below, 
this row structure affects the depinning properties of the 
vortex lattice. 

We next study the transport properties of the differ- 
ent fillings by applying a transport current in either the 
x or y direction. Due to the strong anisotropy of the 
pinning geometry, we find strongly anisotropic depinning 
thresholds in the x direction and f]f in the y direction. 
Figure 0] illustrates the depinning force at the matching 
fields for = /<jx and = fyy. The depinning force in 
the y direction flf is much smaller than f* at all match- 
ing fields, and f% drops nearly to zero above B/B^ = 5 
once interstitial vortices begin to appear in the free chan- 



FIG. 7: Critical depinning force /* at higher fillings B /B$ > 
8. The y axis of the figure has been magnified relative to 
Fig. 4. 



nel region between the pins. At incommensurate fields, 
not shown, we expect the depinning force for both x and 
y directions to be lower than it is at the commensurate 
fields. 

We find a nearly constant f% for B/B^ < 7, while for 
B/B^ > 8, gradually decreases with B/B^. Repre- 
sentative velocity-force curves for x direction driving are 
illustrated in Fig. 03 The curves for B/B^ = 2 through 7 
are essentially identical to the curve shown for B / B^ — 1 , 
with a sharp depinning transition at f% = 22.2 followed 
by an ohmic response. The transition in f£ at B /B^ > 8 
occurs due to a change in the depinning mechanism. For 
B/B^ < 7, vortices initially accumulate inside the pins 
and the pinning troughs upon application of a driving 
force, as illustrated in Fig. EJa) for B/B^ = 7. Depin- 
ning occurs when the pinning force of the pinning troughs 
in the x direction is overcome. When B/B^ > 8, not all 
of the vortices can fit in a single column inside the trough 
because the confining force of the pinning troughs is in- 
sufficient to overcome the vortex- vortex repulsion. In- 
stead, the vortex column buckles in the trough area be- 
tween pins, as illustrated in Fig. EJb) for Bj B^ — 9. De- 
pinning initiates at the buckled points, where the vortex- 
vortex repulsive force is now added to the driving force 
on a buckled vortex; as a result, the depinning threshold 
is considerably lower than for B/B^ < 8 and decreases 
with increasing field. At higher fields the buckling in the 
vortex configurations below depinning shifts from out- 
side the pins to inside the pins, as shown in Fig. [B*^ c) 
for B /B$ = 11, and eventually buckling occurs in both 
locations, as illustrated in Fig.[|Jd) for B/B^ = 17. At 
higher fields, all of the vortices are no longer able to fit 
inside the pinning trough before the depinning threshold 
is reached, and instead the interstitial vortices pile up 
against the repulsive barrier formed by the vortices in- 
side the pinning trough, as shown in Fig.|SJe) and (f) for 
B / B^ = 20 and 27. In this regime, the pinning resembles 
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FIG. 8: Illustration of the vortex motion just after depinning 
under application of a driving force in the y direction (to the 
top in the figure). Black dots indicate vortex positions, pin- 
ning sites are indicated by ellipses, and lines indicate trajec- 
tories followed by the vortices over a period of time. B/B^ = 
(a) 8, (b) 10, (c) 19, (d) 20, (e) 21, and (f) 27. 



extremely strong twinning barriers. 

At B/Bcf, — 9, close to the transition where the ini- 
tial buckling of the vortex configuration before depinning 
shifts from the location between the pins to within the 
pins, we observe a negative dV/dl characteristic as indi- 
cated in Fig. [S] around f% = 16. Here, after depinning the 
vortices initially flow between the rows of pinning sites 
at the location of the buckled points, but as the drive in- 
creases, the flow location shifts and the moving vortices 
instead pass directly through the pinning sites while the 
vortices in the troughs between the pinning sites remain 
pinned. At higher drives there is a step up in V x at 
the point when all of the vortices depin. This feature 
is similar to the dynamic phases observed for symmetric 
periodic pins in Ref. |2(j. For fields above B/B^ = 9, 
depinning initiates at the pinning sites and not in the 
troughs between pins. The initial vortex motion passes 
through the pinning sites, and the vortices in the troughs 
between pins do not depin until higher drives are applied. 
At fields of B/B^ — 17 and higher, when interstitial vor- 
tices are always present even just below the depinning 
transition, the depinning gradually becomes more elastic 
with increasing field, as indicated by the gradual round- 
ing of the velocity force curves at depinning in Fig. 5. 

Depinning in the y direction occurs at much lower ap- 
plied driving forces than depinning in the x direction. 
In addition, when we zoom in on the small driving force 
regime, shown in Fig. Q we observe a nonmonotonic de- 
pendence of fy on field that is caused by differences in 
the vortex configurations. At fillings B/B^ > 7, the 
y direction depinning forces are similar in scale to the 
vortex-vortex interaction forces, permitting the system 
to be sensitive to vortex configuration changes. This is in 
contrast to the case for x direction depinning where the 



Q 




FIG. 9: Finite time total diffusion D t as a function of tem- 
perature for representative vortex fillings B jB$ = (bottom to 
top) 1, 3, 5, 8, 10, 13, 17, and 22. 




FIG. 10: Diffusion in x (D x , bottom curves) and y (D y , top 
curves) directions for fillings BjB$ = (a) 1, (b) 4, (c) 7, (d) 
10, (e) 15, and (f) 27. 



extremely strong pinning forces dominate the response 
of the system. The y direction depinning force f% is en- 
hanced for B/B^ = 10, 20, and 27, as indicated in Fig.[7| 
To understand the origin of this enhancement, we plot 
the vortex trajectories at driving forces just above de- 
pinning for representative fields in Fig. El At B/B^ — 8 
in Fig. Eta)) the vortices in the interstitial channels form 
straight columns and can flow directly in the y direction 
upon depinning. For B/B^ — 10, the buckling of the 
interstitial vortex channels into the spaces between the 
pinning rows disrupts the straight interstitial columns of 
vortices, and the vortex paths upon depinning in the y 
direction wind irregularly. A higher driving force must be 
applied to produce flow through winding channels com- 
pared to the straight channels, since depinning can no 
longer occur via simple propagation of a solitonlike vor- 
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(d) (e) (f) 



FIG. 11: Illustration of vortex motion during the anisotropic melting process. Ellipses indicate the locations of the pins, black 
dots represent vortex positions, and lines show the trajectories of the vortices over a fixed period of time, (a) B / — 4, 
T = 7.05. (b) B/Bj, = 7, T = 3.90. (c) B/B$ = 10, T = 2.46. (d) B/B$ = 15, T = 2.34. (e) B/B^> — 20, T — 1.86. (f) 
B/B<p = 27, T= 2.25. 



tex pulse. 

As the number of vortices inside the interstitial chan- 
nel increases, more than one channel of flow can open, 
as illustrated in Fig. Etc) for B/B^ = 19. This multi- 
ple channel flow occurs relatively easily and is associated 
with a low value of f]!. At B/B^ = 20, where we ob- 
served well formed rows of vortices parallel to the x axis 
in Fig. Etb), the vortex flow is hindered for depinning 
in the y direction since the initial depinning paths wind 
back and forth, as shown in Fig. Et°0 ■ Addition of one 
more vortex per unit cell distorts the rows of vortices 
parallel to the x axis and permits a smoother channel 
flow to occur with lower f%, as illustrated in Fig. EI e ) 
for B/B^ = 21. Finally, we observe an increase in /J 7 at 
the highest field we studied, B/B^ = 27. In this case, 
the columns of vortex motion at depinning do not mix; 
instead, the interstitial vortices move through five paral- 
lel channels, as shown in Fig. EtO- Two of the channels 
contain only a single vortex per unit cell, and these two 
channels help to effectively jam the flow at lower driving 
forces. 

We next consider the melting properties of the vortex 
configurations at zero applied driving force. We measure 
the total diffusion D t = \v(ti) — r(to)\/8t as a function of 
temperature by calculating the distance traveled by the 
vortices during a fixed time period 5t = t\ — to. For a 
fixed choice of 5 t , melting appears as a noticeable change 
in the slope of D t with temperature. Figure ED shows 
representative D t vs T curves for several values of B/B^. 
The melting temperature decreases monotonically with 



B/B+. 

Due to the strong anisotropy of the pinning structure, 
the melting occurs anisotropically. To measure this, we 
compute the diffusion in the x and y directions separately, 
D x = \r x (h) - r x (t )\/5t and D y = \r y (ti) - r y (t )\/5t. 
We show representative plots of D x and D y for six val- 
ues of B/Btp in Fig.Uni At B/B^ = 1, in Fig-HOTal. the 
vortices remain confined within the pins for T < 10, but 
since the pins are more extended in the y direction than in 
the x direction, D y is larger than D x at low temperatures, 
indicating that the vortices undergo larger thermal excur- 
sions within the pins in the y direction. At B/B^, — 4, 
shown in Fig. HOf b). the change of slope in D y occurs at 
lower temperatures than the corresponding slope change 
in D x . In this case, thermal activation causes a small 
number of the vortices to jump out of the wells and oc- 
cupy interstitial sites once T « 6. These vortices remain 
confined within a particular interstitial channel but are 
able to diffuse slowly in the y direction from one inter- 
stitial site to another, as illustrated in Fig. Illf a'). There 
is no significant x direction diffusion until T sa 8, when 
all of the vortices begin to hop in and out of the pinning 
sites. When interstitial vortices are present from the be- 
ginning, such as at B/B^ = 7 shown in Fig. HUT c). D y 
is significantly larger than D x over a sizable temperature 
range. This is indicative of the formation of an intersti- 
tial liquid, where vortices diffuse freely in the y direction 
in one-dimensional channels between the columns of pins, 
but the x direction vortex diffusion is suppressed. A rep- 
resentative example of the interstitial liquid is plotted in 
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FIG. 12: Squares: Melting temperature in the y direction, 
Tfh, as measured from D y . Circles: Melting temperature in 
the x direction, T„, as measured from D x . 
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FIG. 13: Specific heat C v for selected fillings B/B^ = (a) 2, 
(b) 5, (c) 9, (d) 11, (e) 14, (f) 16, (g) 20, (h) 25, and (i) 27. 



Fig. II If b) for B/Bcj, = 7. A plateau begins to form in D x 
for B/B^ = 10 and above, as illustrated in Fig. HOf d). At 
these higher fields, the interstitial channel contains more 
than one column of vortices, and the onset of y direction 
diffusion within the interstitial channel is accompanied 
by a fixed amount of x translation due to the formation 
of ringlike structures in the vortex trajectories, as shown 
in Fig. Hir e) . D x remains fixed at a constant value while 
D y continues to increase with temperature as long as the 
interstitial liquid exists. Similar behavior occurs for two 
interstitial vortex channels at B/B^ = 15, as shown in 
Fig. HTH e) and Fig. UTT d). as well as for three intersti- 
tial vortex channels at B/B^ = 20, plotted in Fig. lllf e). 
and for four interstitial vortex channels at B/ B$ = 27, 
illustrated in Fig.HCTfl and Fig.Egf). 

In Fig. ^] we summarize the behavior of the melt- 
ing temperatures in the x and y directions, T^j and 
T^, as determined from the diffusion. The melting oc- 
curs at high temperatures and is nearly isotropic for 
B I B,p = 1, but anisotropy develops rapidly as the field is 
increased and becomes particularly pronounced once in- 
terstitial vortices appear at B / B$ = 5. The y-direction 
melting temperature, T,^, drops rapidly with field for 
5 < B/Bj, < 10 as the interstitial vortex channel de- 
velops, while T^j changes only slowly with field in this 
regime. Once the interstitial channels have fully formed 
and begin to buckle into a configuration that is two or 
more vortices wide, at B/B^ > 10, levels off and 
begins to drop more rapidly as the space between intersti- 
tial channel vortices and pinned vortices decreases with 
field. Above BjB^ = 15, when the vortex configurations 
inside the pins begin to buckle, also levels off at a low 
value. 

The melting process at higher fillings in general pro- 
ceeds as follows. The first stage is contour melting, in 
which the particles diffuse along the columns of pinning 
and maintain their spacing in columns. Exchange be- 



tween columns happens at points where the columns are 
narrow, which generally falls in between pinning sites. As 
the temperature increases, this channel structure is lost 
and the system forms a series of liquid channels where 
all structure is destroyed. The particles inside the pins 
move freely inside the pins but remain confined. The 
final stage of melting occurs when the pinned particles 
depin and begin to move along the pinning channel at 
the same time as particles begin to jump in the x direc- 
tion from the pinning channel into the unpinned channel 
liquid and vice versa. 

The multiple stages of melting in the two directions 
also appear as signatures in the specific heat, 



C V (T) = ((E(T) 2 ) - (E(T)) 2 )/T 2 



(4) 



obtained from the average energy of the system E(T). 
To measure C v , we begin with aT = vortex configu- 
ration obtained from a slow anneal. We then gradually 
increase the temperature in increments of AT. At each 
temperature we compute the average energy of the sys- 
tem, E, and average E(T) over 20 to 600 realizations 
per filling fraction. We performed the largest number 
of realizations for small values of BjB^. We obtain the 
specific heat C v by taking the derivative of E(T). The 
resulting curves C v (T) are illustrated for selected fillings 
in Fig. El For B/B^ < 5, melting occurs when the 
vortices are thermally activated out of the pinning sites. 
This produces a single step up in C v at relatively high 
temperatures near T ~ 10, as illustrated in Fig. I13f a) 
for B/B^ = 2. Once interstitial vortices appear, as in 
Fig. I13f b~) for B/B^ = 5, a new signature appears in C v 
at a lower temperature T ~ 6, corresponding to the tem- 
perature at which when the interstitial vortices begin to 
diffuse and form an anisotropic interstitial liquid. The 
magnitude of this second peak in C v varies with filling, 
as shown in Fig. 113( d) through (i), and the peak shifts to 
lower temperature as B/B^ increases, consistent with the 
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decreasing melting temperature shown in Fig.^J Due to 
the presence of the pinning, all of the melting transitions 
in this system should be thermally activated crossovers. 

In summary, we used numerical simulations to study 
the configurations, dynamics, and melting properties of 
vortex lattices interacting with elliptical pinning sites for 
both low and high matching fields. The configurations we 
obtain agree well with experimental vortex images from 
the system we are modeling, and as in the experiment 
we find that the saturation vortex number of the indi- 
vidual pins increases with applied field. The depinning 
thresholds are highly anisotropic, particularly at higher 
fields, and the vortex configurations upon depinning are 
very different depending on whether the current is ap- 



plied along or perpendicular to the pinning channels. The 
depinning threshold in the x direction occurs at large 
driving forces and decreases monotonically with B/Bq, 
while the depinning threshold in the y direction occurs 
at low driving forces and shows nonmonotonic behavior 
in response to the vortex configurations at each field. We 
find anisotropic diffusion when the system melts, along 
with a second structure in the specific heat correspond- 
ing to the presence of an anisotropic interstitial vortex 
liquid at intermediate temperatures. 

We thank G. Karapetrov for useful discussions. This 
work was supported by the U.S. Department of Energy 
under Contract No. W-7405-ENG-36. 
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